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We outline a new class of continuous- variable graph states that can be useful to describe entangle- 
ment, and also multimode squeezing, in an optical frequency comb. We show that a particular case 
of such states coincides with the squeezing of the carrier-envelope phase, or of the total intensity, 
of a mode-locked laser. We then discuss the experimental issues related to measuring the quantum 
noise of the carrier-envelope phase of a frequency comb and show that this can be carried out by 
use of quantum heterodyne multiplexing, a technique that may be useful to measure multipartite 
entanglement as well. 
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a. Introduction. The description of multipartite en- 
tanglement by graph states [l| was recently developed 
in conjunction with the one-way quantum computing 
model, in which entangled states of the cluster type are 
used [2]. As the name indicates, a graph quantum state 
is completely described by a mathematical graph, i.e. a 
set of vertices connected by edges [1] as per precise con- 
ventions. A vertex represents a physical system, e.g. a 
qubit (or qubit, 2-dimensional Hilbert space) qudit (d- 
dimensional Hilbert space), or qunat (continuous Hilbert 
space) [3], in a specified quantum state. An edge between 
two vertices represents the physical interaction between 
the corresponding systems. In this article, we will use 
different graph definitions, i.e. conventions for edges and 
vertices, from those used for cluster states. We will also 
use continuous quantum variables, that is, position and 
momentum, or quadrature amplitudes of the quantized 
electromagnetic field (which are equivalent to the former 
by virtue of the quantum harmonic oscillator formalism) . 
Note that continuous- variable (CV) entanglement can al- 
ways be recast in terms of discrete variables in an infinite 
dimensional Hilbert space by use of the Schmidt decom- 
position 0, [Bj], hence qunats can always be viewed as 
qudit s. 

Our group proposed a compact and scalable exper- 
imental implementation of Greenberger-Horne-Zeilinger 
(GHZ) states by use of multimode squeezing generated 
by concurrent optical nonlinear it ies in a single, sophis- 
ticatedly engineered, nonlinear optical medium @, Q. 
A demonstration of such a nonlinear medium was sub- 
sequently made in periodically poled potassium titanyl 
phosphate 0. More recently, we showed theoretically 
that the generation of CV cluster states in a single op- 
tical parametric oscillator (OPO) is possible and gave a 
detailed experimental proposal for a square cluster [9[. 
Our method can be described by simple graph states, as 
explained below. In this paper, we extend the application 
of such graph states to multimode squeezing and then fo- 
cus on the squeezing of the total phase of a set of coherent 
fields, regularly spaced in frequency ("frequency comb"), 
such as emitted by a mode- locked pulsed laser. This to- 
tal phase is also called the carrier-envelope phase and 
its control is crucial in the use of femtosecond-laser fre- 



quency combs for defining frequency standards [TO, EU • 
We then describe experimental procedures for generating 
and for measuring, by quantum heterodyne multiplexing, 
the squeezing of the carrier-envelope phase. The latter 
being a CV GHZ entanglement witness, these procedures 
can be used for entanglement measurements as well. 

b. Two-mode-squeezed graph states. The generation 
of scalable CV GHZ states was first proposed by van 
Loock and Braunstei n 11 2| and demonstrated experimen- 
tally for 3 modes [l|,[l|. In this method, N single-mode 
squeezed states are generated by N independent optical 
parametric amplifiers (OPA) and subsequently interfered 
at a 2A/"-port beam splitter network with all optical paths 
controlled. A similar setup can be used to generate linear 
[IH or even arbitrary [16] cluster states. Such proposals 
nonetheless lead to nontrivial, if polynomial, experimen- 
tal complexity when scaling to large N. In contrast, our 
approach uses a single OPA, or OPO, as a multimode 
squeezer that delivers N entangled modes directly, pro- 
vided the nonlinear medium constituting the OPA phase- 
match simultaneously all two-mode interactions between 
every one of the optical modes Q . This lends itself quite 
naturally to a graph description. As an example, we con- 
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FIG. 1: GHZ graph states with two-mode squeezing interac- 
tions only. The N modes are denoted by the vertices, con- 
nected by N(N - l)/2 interaction edges, (a) N = 2, (b) 
N = 3, (c) N = 4, (d) N = 5. 

sider GHZ states depicted in Fig.[TJ These are created by 
the second-order nonlinear interaction Hamiltonian [6] 

H t =ihn^a ] n al + H.c. (1) 

where k is the nonlinear coupling constant and each term 
is represented by an edge of the graph. Each mode is 
a vertex and a qudit in the photon-number basis or a 



qunat in the quadrature amplitude basis labeled by 
frequency as well as the wave vector and polarization di- 
rections. The set of interacting modes is defined by the 
resonant modes of the OPO optical cavity, inside which 
the nonlinear concurrent emitter is placed. It is straight- 
forward to see from the Heisenberg equations for the 
modes {a n }i...jv that the matrix describing this system of 
coupled differential equations is proportional to the ad- 
jacency matrix [1] of the graph so defined. In this case of 
a GHZ state, the graph is fully connected and the adja- 
cency matrix A has elements A k \ = l — 5ki, where S is the 
Kronecker symbol, has eigenvalues {TV — 1,— 1,...,— 1} 
and respective (unnormalized) eigenvectors {(1,1,...,1), 
(1,-1,0,...,0), (1,0,-1,0,. ..,0), (1,0,. ..,0,-1)} 
[(| . In terms of multimode squeezing of the field quadra- 
tures Xj = ctj + ctj ("amplitude") and Pj = z(aj — ctj) 
("phase"), this eigensystem yields the following solution 
to the quantum evolution: the phase sum P = J2f=i Pj 
is squeezed P(t) = P(0)e _( ^ Ar_1 ^ £ (note the squeezing 
enhancement factor N — 1), along with all amplitude 
differences X kl = X k - X u X kl (t) = e~ Kt X kl (0) p^. 
In the limit of infinite squeezing, this means the result 
of the quantum evolution is a common eigenstate of all 
these operators, which coincide with the entanglement 
witnesses of a GHZ state [HI, Qil . Such a state has the 
unnormalized expressions, in the amplitude-quadrature 
basis, 

J \x) 1 \x) 2 • • • \x) N dx. (2) 

We have named this particular topology of nonlinear in- 
teractions a nonlinear optical concurrence (NOC) — not 
to be mistaken for the concurrence as a mathematical 
measure of entanglement [20| — because it mandates the 
co-occurrence of several nonlinear interactions (edges) for 
the same mode (vertex). The N = 2 entangling NOC has 
a trivial graph and is already well-known from previous 
studies of the Einstein-Podolsky-Rosen (EPR) paradox: 
the corresponding entangled EPR state can, indeed, be 
generated from either a single nondegenerate OP A 12 1| o r 
from the interference of two single- mode OPAs [22|, [23[ • 
Implementing an entangling NOC would seem a fairly 
involved task, for reasons we detailed previously @, 0], 
but we have proven it feasible [8]. One issue is that CV 
GHZ graphs cannot have any self-loops of the same type 
of interaction lest entanglement be lost. Indeed, as we 
previously established [6] , the protocol of van Loock and 
Braunstein [12[ corresponds to the Hamiltonian 

h 2 = ih K (j2 a H - ^ E °r) + H - c - (3) 

\n<p n J 

which corresponds to a self-loop graph because of the 
second, degenerate interaction term. However, the graph 
is now weighted: all self- loops describe an interaction of 
opposite sign and different strength (—(N — 2)/2) with 
respect to the 2-vertex edges (Fig. [2]). Experimentally, it 
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FIG. 2: GHZ graph states with single-mode (degenerate, 
loops) as well as two-mode (nondegenerate, 2-vertex edges) 
squeezing interactions. Here the loops all correspond to in- 
teractions opposite to those of the edges, i.e. downconverting 
versus upconverting, or vice versa. 

means that, if edges are parametric downconverting in- 
teractions, then the self-loops must be upconverting ones, 
or vice versa. Such a configuration is clearly incompati- 
ble with a single OPO, which is why the NOC scheme of 
Fig. [1] is interesting for the single-OPO implementation. 

If one now changes the signs of the self- loops, however, 
the Hamiltonian becomes 

H 3 = ih K (j2 a H + E t~ ) + H - c -' ( 4 ) 

\n<p n J 

with the corresponding graphs given in Fig. [3l The ad- 
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FIG. 3: Graph states for total-phase or total- amplitude 
squeezing. 

jacency matrix has all its elements equal to 1 and has 
eigenvalues {N, 0, • • • , 0} with the same respective eigen- 
vectors. The only squeezed operator here is the phase 
sum P(t) = P(0)e~ NKt (note the squeezing enhancement 
factor N). The other entanglement witnesses, the ampli- 
tude differences, are now constants of the motion (zero 
eigenvalue), that is, not squeezed and hence no longer 
diagonalized by the resulting state in the limit of infinite 
squeezing. 

c. A phase-sum squeezer: the broadband degenerate 
OPO. The Hamiltonian H3 is actually very easy to im- 
plement, compared to H\$'- it describes an OPA or OPO 
pumped by a train of mode-locked pump pulses whose 
repetition rate is equal to the free spectral range of the 
OPA cavity, as sketched in Fig. [U It is then easy to 
see that each pump frequency in Fig. 0] contributes one 
skew diagonal of ones to the adjacency matrix, thereby 
yielding the aforementioned uniform adjacency matrix, 
corresponding to the Heisenberg system for H%. Note 
that the match between the pump's repetition rate and 
the OPO cavity free spectral range is crucially impor- 
tant but can be easily controlled by use of an electronic 
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FIG. 4: Experimental implementation of the graph states 
such as Fig. [3] All pump field spectral components have the 
same polarization. The open and filled circles label, respec- 
tively, the pump components who have a degenerate inter- 
action with the signals, or not. For clarity, not all signal 
couplings have been drawn. 



servo loop. For pump and signal spectral ranges within 
the phase-matching bandwidth of the OPO crystal, one 
should therefore obtain phase-sum squeezing of the signal 
comb. (Note that the OPO cavity linewidth is typically 
much narrower than the phase-matching bandwidth and 
than the free spectral range.) 

In the case of a broad femtosecond frequency comb 
whose high-frequency end of the spectrum would be the 
pump modes, and whose low- frequency end of the spec- 
trum would be the signal modes, it is then not clear 
whether one would obtain carrier-envelope phase squeez- 
ing since only the sum of the signal phases would be 
squeezed, according to the analysis of this paper. How- 
ever, one must keep in mind that this paper assumes 
undepleted classical pump modes, which is not necessar- 
ily true in that situation. One should then include the 
quantum behavior of the pump fields during the interac- 
tion, using more sophisticated quantum optical methods 
[24|. Such a study is compelling but out of the scope of 
the present paper. Note that a complementary mathe- 
matical analysis of the same system was carried out in- 
dependently of this work [25] but the squeezed physical 
observables were called "supermodes" and not identified. 

d. Phase-sum squeezing detection: quantum hetero- 
dyne multiplexing. Considering again the case of inde- 
pendent pump and signal combs, we now show that the 
detection of the signal phase-sum can be achieved rela- 
tively easily by balanced heterodyne detection. Hetero- 
dyne detection does not usually allow squeezing to be 
measured because it lets in vacuum fluctuations. How- 
ever, there exists a particular case where this effect can 
be occulted. This case is the balanced heterodyne detec- 
tion of two frequency nondegenerate modes by use of a 
single local oscillator (LO) field of frequency the exact 
average of the modes' frequencies. This is the method 
that was used in the first squeezing measurement (26| . 
In the photocurrent difference signal, all DC signals are 
filtered out, as well as the interference signals between 



modes that impinge onto the beam splitter through the 
same input port. We consider the configuration where all 
quantum squeezed modes are input into the same beam- 
splitter port and the local oscillator is sent into the other 
port, as depicted in Fig. [5l 




FIG. 5: Detection of total-phase squeezing. The a input 
port receives the 2 modes whose phase-sum is squeezed. The 
b mode receives a coherent LO field. The vacuum modes 
that contribute to the signal at frequency Q are identified by 
dashed lines. The 2 vacuum modes at the extremities of the 
b input spectrum are the "image" modes. 

If we denote by aj the annihilation operator of the a 
input mode at frequency uo + jO, and bj that of the b 
input mode at frequency uj + jO, the operator detected 
by the setup of Fig. [5] after photocurrent subtraction is 

S = a)b + atf = E a JE^ + E a iE 6 I ( 5 ) 

j k j k 

~ (aLi + a\) b + (a_i + ai) 6q (6) 



(1+) 



(7) 



in the limit of a classical LO whose amplitude is much 
larger than those of the squeezed fields, in which case 
the beats of the a fields with the image bands in b can be 
neglected. We have posed bo = ftoe %e a complex number, 
the quadrature A 



-id, 



e ie a], and a^'+) = a. 



-j ' 



For the LO phase 6 = 7r/2, one will measure, at 
photocurrent frequency the sum of the signal phases 

It is straightforward to extend the original two-mode 
measurement to a larger number of modes by making 
use of the frequency degree of freedom of the photocur- 
rent: in Fig. [6l each mode pair a±j of the signal comb 
is beat against the LO, yielding a signal at (2j — l)fl in 
the photocurrent difference. By virtue of the properties 
of balanced heterodyne detection stated above, demodu- 
lating the photocurrent difference at all such frequencies 
and summing all results will yield the operator 



S = f3 J2A j+) =f3 £ A<j\ 
j'=i 



(8) 



j =—n 



which coincides with the total phase P for = tt/2. In 
the expression above, note that n ideally equals the total 
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FIG. 6: Detection of total-phase squeezing for a frequency 
comb. 

mode number (TV) in the squeezed operator. In practice, 
n <C N is possible due to a detection bandwidth limit 
of nQ. This would create a problem since the fact that 

Y^^L-n Pj i s squeezed does not imply that Y^j=- n Pj * s 
squeezed. This problem can easily be surmounted, how- 



ever, by the use of a comb of several LO fields, spaced in 
frequency by nVt (assuming the detection band-shape is 
a rectangular window). The total output signal can then 
be demodulated and summed in exactly the same man- 
ner to yield the same result. Overall, this use of the fre- 
quency domain for encoding quantum information more 
densely is a quantum analog of FM radio. We name this 
detection method "quantum heterodyne multiplexing." 

In conclusion, we have shown that the topology 
of graph states, denned with squeezing rather than 
controlled-phase interactions connects to entangle- 
ment and multimode squeezing, and in particular phase- 
sum squeezing, of interest to improve the metrological 
performance of frequency combs. Further theoretical 
study is needed at this point to investigate the possibili- 
ties for the self-referenced metrological comb. Finally, we 
proposed a compact multimode detection method for the 
total phase, which can be useful for measurements of mul- 
tipartite entanglement as well. This work was supported 
by NSF Grants No. PHY-0555522 and No. CCF-0622100. 
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